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ABSTRACT 

Based on results from cold dark matter N-body simulations we develop a dynamical 
model for the evolution of subhaloes within group sized host haloes. Only subhaloes 
more massive than 5 x 10 s M Q are considered, because they are massive enough to 
possibly host luminous galaxies. On their orbits within a growing host potential the 
subhaloes are subject to tidal stripping and dynamical friction. At the present time 
(z = 0) all model hosts have equal mass {M vir — 3.9 x 10 13 M Q ) but different con- 
centrations associated with different formation times. We investigate the variation of 
subhaloe (or satellite galaxy) velocity dispersion with host concentration and/or for- 
mation time. In agreement with the Jeans equation the velocity dispersion of subhaloes 
increases with the host concentration. Between concentrations ~ 5 and ~ 20 the sub- 
haloe velocity dispersions increase by factor of ~ 1.25. By applying a simplified tidal 
disruption criterion, i.e. rejection of all subhaloes with a tidal truncation radius below 
3 kpc at z = 0, the central velocity dispersion of the 'surviving' subhaloe sample in- 
creases substantially for all concentrations. The enhanced central velocity dispersions 
in the surviving subhaloe samples are caused by a lack slow tangential motions. Ad- 
ditionally, we present a fitting formula for the anisotropy parameter which does not 
depend on concentration if the group-centric distances are scaled by r s , the character- 
istic radius of the NFW-profile. Since the expected loss of subhaloes and galaxies due 
to tidal disruption increases the velocity dispersion of surviving galaxies, the observed 
galaxy velocity dispersion can substantially overestimate the virial mass. 
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1 INTRODUCTION 

Clusters of galaxies have been extensively studied in all 
available wave bands. They belong to the most mas- 
sive and eye-catching structures in the Universe and are 
thought to closely resemble the cosmic composite of 85 
per cent dark mat ter and 15 per cent baryonic matter 
ijSpergel et alj |2006h According to the hierarchical model 
for structure assembly, these systems preferentially have 
formed late. For dark matter haloes it has been shown 
that the formation redshift i s correlated with con centra- 
tion (e.g. iBullock et all l200ll : lw"echsler et all I2002T ) there- 
fore clusters present a class of low concentration objects. 
In contrast, groups of galaxies span a whole range of for- 
mation times and concentrations. Some resembl e the prop- 
erties of clusters (IZabludoff fc Mulchaevl Il998l : iMulchaevI 
l200d : IZabludoff fc Mulchaevl l200Ch . having low concentra- 
tions, others show an opposite behaviour. 

Also the properties of the galaxy populations vary 
strongly amongst individual groups. Some groups have lumi- 



nosity functions similar to those observed in clusters but 'fos- 
sil groups' display a central giant elliptical g alaxy with only 
few intermediate-lumi nosity galaxies (e.g. iPonman et al.l 
Il994l :lj ones et al]|2003l ). The optical properties of the central 
bright galaxies are consistent with an origin as the merged 
combination of a more typical, less over- luminous, galax- 
ies. The variations amongst the satellite populations are 
thought to originate from varying fo r matio n times. Using 
N-body experiments iD'Onghia et al.l (|2005f ) demonstrated 
that the earliest formed group-sized dark matter haloes 
reve al structural similarity to observed fossil groups (see 
also ISommerTl arscri 200g). The early formation times of 
fossil groups also impacts satellite galaxies that do not 
merge with the central giant. Based on the observed deple- 
tion of intermediate -luminosity galaxies in the fossil group 
RX J1552. 2+2013, iMendes de Oliveira et al] (|2006l 1 sug- 
gests that the me mbers of this system may be a f fected 
by tidal disruption. |Pa Rocha fc Mendes de Oliveiral (120051 ) 
find substantial amount of intra-group light in two of three 
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compact groups which may, at least partly, come from tidally 
dissolved satellite galaxies. These observations indicate a 
substantial variation of the satellite population with time 
either by merging or by tidal dissolution. 

In addition to the concentration-age and the luminosity- 
age relations, we describe here a concentration-velocity re- 
lation. According to the Jeans equation, galaxies orbiting in 
more concentrated host haloes have higher velocity disper- 
sions. This effect is most pronounced in the central regions 
of the host haloe. As discussed above, early formed groups 
are more concentrated, leading to the conclusion that the 
central satellite galaxies in fossil groups move faster than 
central galaxies in less evolved groups of comparable size. 
Additionally, if tidal disruption reduces the number of satel- 
lite galaxies, the velocity dispersion of the remaining galax- 
ies is expected to increase even higher. The last statement 
is supported by t he following considera tion. Analysing N- 
body simulations, iDiemand et al.l l|2004h find an enhanced 
velocity dispersion of subhaloes compared to the disper- 
sion of the diffuse dark matter component. Subsequently, 
iFaltenbacher fc Diemandl i|200rj) argued that the loss of a 
slow moving, earlier accreted, subset of the subhaloe pop- 
ulation is responsible for this bias. In analogy, if satellite 
galaxies are prone to tidal dissolution, then early accreted, 
slow moving, galaxies are preferentially dissolved. The lack 
of low velocity galaxies causes an enhanced velocity disper- 
sion of the remaining galaxies. An observational hint for such 
a scenario may be the rec ent observation of the fo ssil sys- 
tem RX J1416.4+2315 by l|Khosroshahi et alJlioOrj) . Based 
on 18 member galaxies they find a velocity dispersion that 
is nearly twice as high as expected from the X-ray analysis 
of this system. 

These findings motivated us to study the relation be- 
tween the host concentration and dynamics of the subhaloe 
populations. The paper is organised as follows. § [2] intro- 
duces the dynamical model used to trace the evolution of 
substructure in host haloes of different mass accretion his- 
tories. § [3] deals with substructure as point-like particles, a 
reference for the subsequent more complex investigation. §[3] 
presents dynamical results when substructures are treated 
as extended units prone to tidal stripping and dynamical 
friction. Finally, § concludes with a summary. 



2 DYNAMICAL MODEL 

We aim to investigate the correlation between the host con- 
centration and subhaloe dynamics. In principle cosmologi- 
cal N-body simulations can best tackle these kinds of prob- 
lems. However, such simulations should cover a large volume 
for satisfactory host haloe statistics and, at the same time, 
provide high resolution for the accurate determination of 
the subhaloe dynamics. Here we pursue an alternative ap- 
proach. We employ semi-analytical models to gain a basic 
understanding of the processes which shape the subhaloe 
distributions within equal mass hosts of various concentra- 
tions. 

In recent years high resolution CDM simulations have 
substantially improved our understanding of cosmic struc- 
ture formation. This insight can be used to model the dy- 
namical evolu tion of subha l oe in a semi-analytical man- 
ner. Recently, IZentner et al.l (|2005l ) demonstrated that this 



kind of semi-analytical modelling achieves good agree- 
ment with state of the art N-body simulations. Simi- 
lar approache s have been used for a variety of scientific 
goals (see e.g. iBullock et al.ll2000l.l200ll; IZentner fc Bullock 
20031; iKoushiappas et al.ll2004l; iTavlor fc Babulll200ll . 120041 

Islam et al.ll2003l ; Ivan den Bosch et al.ll2005l). 

Here we extend the models present e d in | Mathews et alj 

l|2004h and IFaltenbacher fe Mathewsl (|2005f) . Originally, 
these models were designed to reproduce the well observed 
number density distribution of satellite galaxies of the NGC 
5044 group. These models have been modified for the pur- 
poses of present investigation. Basically, the model inte- 
grates the orbits of subhaloes from the accretion onto the 
host to the present (z — 0) within the deepening gravita- 
tional potential of the host haloe. By means of this the ori- 
gin of the present subhaloe distribution can be investigated. 
Also, the code allows us to assign a mass profile to the en- 
tering subhaloes which in turn can be used to compute the 
modification of their orbits by tidal mass loss and dynamical 
friction. 



2.1 Host haloe 

Based on the analysis of N-body simulations Wcchslcr et al. 



found that the mass growth of CDM haloes can be 
described by 



M v (a) = M vfi e- 2a f(i-^ 



(1) 



where a = l/(l + z) and a,f = 1/(1 + Zf ) is the cosmic expan- 
sion factor at the formation redshift 2/. Given the present 
virial mass of a haloe, the mass accretion history is fully de- 
termined by the only remaining free parameter a/. Individ- 
ual accretion histories may differ strongly from this descrip- 
tion, but on average the mass accretion is well approximated 
by this formula. Most importantly, the mass built-up of a 
haloe is separated into two distinct phases, namely an early 
rapid accretion phase and a subsequent period of modest 
accretion. Knowing M v (a) also allows to compute the evo- 
lution of the virial radius R V i r (a). According to the spherical 
collapse model the virial radius R v i r (a) is the radius which 
includes a mean density of A(a)p c (a), where A(a) is the 
viri al over-density in units of the critical de nsity p c (a) (see 
e.g.lEke" et al. 1998|; iBrvan fc Normanlll998l ) . 



Rvir 



/ 3M(o) 



1/3 



(2) 



V 4tt A(a)p c (o) 

Anticipating an NFW density profile l|Navarro et al.lll997f ). 
the host mass and concentration are s ufficient to completely 
determine its properties. According to lWechsler et al"] (|2002l 'l 
the evolution of the concentration is given by 



Cvir — 



(3) 



where ci = 5.125 according to IZentner et al.l (|2005T ) and 
ao, the cosmic expansion factor at the time of observation, 
equals 1 in this context. In our model Eqs. Q] and [3] com- 
pletely characterise the growth of the host and its gravita- 
tional potential which in turn determines the orbits of the 
subhaloes. At z — all host haloes approach the same virial 
mass (M Vi q = 3.9 x 10 13 Mq), but they exhibit concentra- 
tions in the range between 2 and 20 (integer steps, resulting 
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in 19 output sets). The low concentration haloes are cur- 
rently in the rapid accretion phase. According to Eq. [3] the 
formation of the least concentrated haloes takes place in fu- 
ture times. The most concentrated haloes formed at redshifts 
z ~ 3. The evolution of the host haloe as described here 
does not explicitly include merging events. Such events may 
rearrange the phase space distributions of the subhaloes. 
However, the rapid accretion phase at early times causes 
the fastest change of the potential over the whole accretion 
history of a haloe. To a certain extent these violent merg- 
ing processes are covered by the steep slope of the accretion 
(equation [TJ at Zf. 

In summary, the mass growth of the host haloes mod- 
elled according to Eq. [T] provides the background potential 
for the dynamical evolution of the subhaloes. All hosts have 
a mass of M„,o = 3.9 x 10 13 Mq at z — and concentrations 
ranging from 2 to 20 (integer steps). According to Eq(3] con- 
centrations are associated with formation times, resulting in 
different accretion histories for the host haloes despite their 
equal masses. 

2.2 Subhaloes 

After introducing the model for the host we now descirbe 
the host's substructure population. Two samples of 100000 
orbits are traced for every host haloe. The first sample, re- 
ferred to as the point mass sample, follows the evolution of 
point-like particles within the potential of the host. In the 
second sample, dubbed the subhaloe sample, the orbits are 
modified by the effects of tidal stripping and dynamical fric- 
tion which involve the attribution of density profiles to each 
substructure. The point mass sample is mainly used to test 
the agreement of the model with N-body simulations on a 
particle basis. The subhaloes in the subhaloe sample have 
masses > 5 x 10 8 Mq which renders them possible candidates 
for hosting galaxies. Obviously, the sum of all 100000 sub- 
haloe masses exceeds the mass of the host haloe. Therefore, 
the results for the subhaloe distributions should be consid- 
ered as the outcome of stacking many identical groups. In 
the next two paragraphs we describe the generation of the 
initial conditions and initial properties of the subhaloes in 
the subhaloe sample. 

2.2.1 Orbital initial conditions 

We assume that the initial space density of the subhaloes 
is proportional to the dark matter, i.e. subhaloes enter the 
host haloe in proportion to the dark matter. To generate 
this particular accretion pattern the individual arrival times 
of the subhaloes are assigned by a random process using the 
mass of the host as a proxy for the cosmic time. 

Mi = RiM„, (4) 

Here < TZa < 1 is a random number and M„,o the virial 
mass of the host at z = 0. By setting Mi = M (a) (Eq. [TJ 
Mi can be converted into the cosmic expansion factor a,i or 
equivalently zt, which is considered as the accretion time of 
the ith subhaloe. Subsequently Eq.[2]is used to derive virial 
radius of the host, r», at that time. The ith subhaloe is as- 
sumed to be gravitationally accelerated by Mi from stand- 
still at the turnaround radius (r tU m = 2rj) until it reaches 



the virial radius, n. The resulting absolute values of the 
velocity at n is 

= (GMt/r t ) 1/a (5) 

which also equals the circular velocity at fj. To assign an 
angular momentum to each subhaloe we invoke another ran- 
dom process designed to reproduce the distribution of circu- 
larities found in N-body simulations. The circularity is de- 
fined as e = J I J c where J is the angular momentum of the 
orbit and and J c is the angular moment um of a circular orbit 
with the sam e orbital energy (see e.g. lLacev fc Cole]|l993l ; 
iTormenll 19971 ). IZentner et alj (|2005T l present the orbital cir- 
cularity distribution for subhaloes in N-body simulations at 
the time of entering the virial radius of the host. They also 
provide a fitting formula for the distribution and state that 
it is independent of accretion redshift and subhaloe mass. 
Here we model the probability distribution of circularities 
with a simple quadratic, dp/de = 6e(l — e), w hich is in fairly 
good agreement with the data presented by IZentner et al.l 
(2005) and has a mean circularity of (e) = 0.5. The initial 
angular momentum of a subhaloe is randomly drawn from 
this distribution, 

K L = [ 6e(l - e)de = 3e 2 - 2e 3 , (6) 
Jo 

where TZl is a random number between zero and one. Note, 
the distributions of circularities were derived for subhaloes 
in N-body simulations, here we implicitly assume that the 
same regularity holds for the point mass sample as well. 

Knowing the evolution of the host potential, the accre- 
tion time and the initial orbital parameters of each subhaloe, 
the integration of the orbits can be pursued. Details about 
the integration are presented at the end of the following 
paragraph. Two random processes are involved to generate 
the orbital initial conditions, namely TZa for the arrival time 
and TZl for the initial angular momentum. We assume, that 
the orbital energy is uniquely determined by the host virial 
mass and radius at the time of arrival, which may not be 
exactly the case (see lZentner et al"1l2005l ). However, for sim- 
plicity we do not implement an additional mechanism to 
spread the distribution of initial orbital energies. 

2.2.2 Subhaloe properties 

As long as the subhaloes are treated as point-like particles 
and dynamical friction is ignored, the orbits are indepen- 
dent of the subhaloe masses. For a more realistic model we 
must consider the spatial extension of the subhaloes and 
as a consequence the susceptibility to tidal truncation. For 
that purpose a mass is attributed to every arriving sub- 
haloe by a random process which is set up to to generate a 
power law mass function as observed for the field haloes 
in N-body simulations. We adopt dMIdN = M" wher e 
a = -1.86 l|Somerville fc Primackl 119991 : iReed et all 120031 ). 
Consequently, the mass of the ith haloe rm at arrival is given 
by 

m, = [Mf - TZ M (M? -M%)f' x (7) 

where Rm is a random number between and 1, x = 1 + a. 
The lower mass limit is Mi = 5 x 10 8 Mq and the upper mass 
limit is M2 = Mi, where Mi is the current host mass (see 
Eq(4]| . We only focus on the evolution of dark matter haloes 



4 A. Faltenbacher and W.G. Mathews 



above M\ — 5 x 10 s Mq because they are massive enough to 
potentially host visible galaxies. The upper mass limit is cho- 
sen to avoid merging of haloes more massive than the current 
host. Note, since the mass function is so steep our results do 
not depend on the exact choice of the upper mass limit. This 
restriction assures that we follow the host merging tree along 
the most massive progenitors. Ignoring the possible devia- 
tions caused by the centrally conden sed baryonic matter , 
we associate an NFW density profile (|Navarro et al.lll997l ) 
to every arriving subhaloe. We assume that the formation 
redshift of all (sub)haloes is z au bj = 7 (which corresponds to 
daubj = 0.125). By setting the arrival time eti = ao in Eq.0 
the concentration Cj of the subhaloe is computed. The den- 
sity profile of the subhaloe is completely determined by m; 
and Ci. We assume that on the subseque nt orbits the shape 
of th e density profile does not change fsee lKazantzidis et al.l 
|2004 ) however it can be truncated due to tidal forces. The 
tidal truncation radius is estimated by the Jac obi limit (see 
iBinnev fc Tremaind[l987l : lHavashi et all 120041 ), 



r j = 



( m c ,i 

\m(< d) 



1/3 



D 



(8) 



Here D is the instantaneous distance of the subhaloe from 
the centre of the host, M(< D) is the host mass within D 
and m c ,i is the current mass of the subhaloe. If r j is smaller 
than the current radius of the subhaloe, all the matter out- 
side of r j is assumed to be stripped instantaneously and the 
current haloe radius is replaced within rj. 

The orbits of the subhaloes are found by solving 



dr 

dt 



and 



du 



GM(r) r 



du\ 



(9) 



where M(r) is the NFW mass within radius r of the host. 
Orbits are computed from redshift Zi (redshift at arrival) to 
the present epoch, z = 0. For the integration we use a fourth 
order Runge-Kutta integration scheme with adaptive time 

stepping. 

The deceleration b y dynamical friction (|Chandrasekharl 
ll943l ; IColpi et al.lll999h is described by 

— ) = -u47rlnAG 2 m / 9it" 3 [erf(X) - -A-Ae^'] . 

dtj df 7T 1 / 2 

(10) 

Here p(r,t) is the local density of the host haloe, u = u 
is the velocity of the subhaloe, X — m/(V2<x) and cr(r, t) is 
the mean velocity dispersion of diffuse dark matter in the 
host halo. We assume that the total dispersion can be ap- 
pr oximated by the co ld dark matter dispersion pres ented 
by iHoeft et all l|2004h (see also iMathews et~al1 12004| ). We 
choose In A = 3 as suggested by Zhang et al. ( 20021 ) which 
is the point mass approximation. This approximation is jus- 
tified by the very efficient tidal mass stripping which rapidly 
truncates the radius of the subhaloe to a few per cent of the 
virial radius of the host. 

In summary, a mass is randomly assigned to each ar- 
riving subhaloe. The random process is designed to recover 
the mass functions of field haloes found in N-body simula- 
tions. Additionally a concentration is assigned to the sub- 
haloe assuming a common (sub)haloe formation redshift of 
z = 7. Based on mass and concentration an NFW density 
profiles is associated with each subhaloe. The profile is as- 
sumed to maintain its shape on the subsequent orbit within 



the host potential well, however tidal truncation is allowed 
to take place. The subhaloe gets most severely truncated 
at its peri-centre passage. The remaining mass remains con- 
stant until the subhaloe approaches the centre for the next 
time. Between two subsequent peri-centre passages the sub- 
haloe loses some orbital energy due to dynamical friction, 
the peri-centre distance shrinks with each passage, reducing 
the truncation radius. Of course, lower masses experience 
less dynamical friction. The mass of the subhaloes decreases 
in a step-like fashion, as also observed in N-body simul ations 
fe.g lKravtsov et al.|[20o3 : iBovlan-Kolchin fc Mall2006l ). 

In order to investigate how tidal reduction affects the 
properties of the subhaloe population we will construct 
a sub-subhaloe sample containing only subhaloes having 
present truncation radii above the rejection radius r re j > 
3kpc. This choice is somewhat arbitrary, but i n agree ment 
with values given in iFaltenbacher fc Mathewsl ([2005) who 
achieve good agreement between the outcome of a very sim- 
ilar dynamical model for the observed galaxy group NGC 
5044. Subsequently, we refer to the two subsamples as then 
complete and reduced samples, emphasising that the latter is 
a subset of the former. Note, disappearance when r < r re j 
in the reduced subhaloe sample does not necessarily mean 
total physical disassembly of the object, it may just mean 
that it falls below a given detection limit. Such limits are 
apparent in both simulations and observations. The detec- 
tion of subhaloes in simulations is constricted by a minimum 
particle number limit for substructures and the detection of 
satellite galaxies in groups depends on the sensitivity of the 
observations. 



3 RESULTS FOR THE POINT MASS SAMPLE 

The point mass sample follows the evolution of point-like 
particles within the growing potential well of equal-mass 
host haloes with varying accretion histories. Mass depen- 
dent mechanisms, like tides and dynamical friction, are ex- 
cluded. The mass of the particles is irrelevant. In this section 
we (1) demonstrate the agreement of the model with estab- 
lished results from N-body simulations and (2) discuss the 
concentration-sigma relation for this simplified approach as 
a reference for the results for the (complete and reduced) 
subhaloe samples presented in the following section. 



3.1 Accretion histories 

The particles, or later subhaloes, enter the host haloe in pro- 
portion to the dark matter. The growth of the dark matter 
host haloe is described by Eq.[T] The generation of the initial 
condition for each particle includes two random selections, 
the arrival time and the initial angular momentum. After 
the first penetration of the host virial radius the orbits are 
integrated numerically. Fig. [T] displays the expansion factor 
at accretion time versus the current particle distance from 
the centre the host, showing distributions for host concen- 
trations of 6, 10 and 20. Subsequently, whenever properties 
of individual hosts shall be exemplified we use these tree 
concentrations and refer to the accordingly, as the c§, cio 
and C20 hosts. Despite the angular momentum distribution 
imposed on the particles, caustic-like features are apparent 
as discussed in the context of the secondary infall model by 
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Figure 1. Scale factor at time of accretion versus current sub- 
haloe distance to the centre for the three example host haloe 
concentrations of 6, 10 and 20 at z = 0. For clarity, only 5 per 
cent of the 100000 particles, randomly selected from the complete 
samples, are displayed. Even if weakened due to the imposed an- 
gular momentum distribution, there is a clear indication for the 
appearance of caustics as described in secondary infall model by 
Bertschinger (1985). The caustic of the second turnaround is lo- 
cated at or beyond the current virial radius. The caustic of the 
most concentrated host extends furthermost, indicating the high- 
est orbital energies. 




expansion factor 



Figure 2. Accretion histories for three host haloe concentrations 
at z = compared to the analytical description in Wechsler & al. 
(2002), thin solid lines. The negative deviation from the Wechsler 
formula at late times results from the loss of particles in of the 
secondary caustic that extends beyond the current virial radius, 
see Fig. [T] 



iBertschingen (|l985l ). The most prominent feature is caused 
by the caustic of the second turnaround generated by the 
particles when they reach the apo-centre after their first 
peri-centre passage. Invest i gating the dynamics of the NGC 
5846 group lMahdavi et alj (|2005h find observational support 
for the appearance of the caustic of the second turnaround. 
The location of the second caustic is correlated with the 
host concentration. For the ce host the caustic is located 
at the virial radius whereas in the C20 host it is found at 
~ 1.3 times the virial radius. For all host concentrations, 
the second turnaround is experienced by particles accreted 
at a ~ 0.8. 

Fig- [U compares the mass growth of the host haloe (thin 
solid lines) to the particle accretion history for all particles 
located within the host virial radius at z — (three heavier 
lines). Both quantities, host mass and particle numbers, are 
normalised to unity. The deviations at late times are caused 
by particles with orbits currently beyond the virial radius 
also seen in Fig. [T] 

The positive correlation between the group-centric dis- 
tances of the caustics and the concentrations in Fig. [T] arises 
from the more rapid mass assembly of high concentration 
hosts. The thin solid lines of Fig. [2]indicate that at a = 0.8. 
The C20 host has already assembled 90 per cent of the final 
mass compared to only 60 per cent for the c§ host. Thus 
at the time (a ~ 0.8) when the particles are accreted that 
presently pass through the second turnaround the C20 host is 
a factor of 1.5 more massive than the c& host. Therefore its 
virial radius is larger by a factor of 1.5 1 ' 3 = 1.14. However, 
Fig.[T]shows a factor of ~ 1.3 difference between the location 
of the C20 and the eg turnaround. This additional shrinking 
of the second turnaround radius within the less concentrated 
host is caused by the rapid mass growth of the ce host at 
the present time. Thus particles moving outward after their 
first peri-centre passage feel a much deeper potential well, 
which prevents them from reaching their original starting 
point, the first turnaround when they decouple from cosmic 
expansion. In contrast, the highly concentrated hosts don't 
show such a rapid deepening of its potential at late times, 
so the particles can move further out. 



3.2 Density profiles 

The host mass grows according to Eq. [T] and the evolution 
of the concentration is given by Eq. [3] The rate that parti- 
cles (the same applies to subhaloes) enter the virial radius 
of the host is given by Eq. [7] About half of all the particles 
are accreted onto the host before z — 0.6. Of course this 
value varies depending on the accretion histories which in 
our model are determined by the current concentration of 
the host. Fig. [3] compares the number density profiles for 
the three concentrations (c = 6, 10, 20) with NFW-profiles. 
The agreement of the analytical profiles with the results of 
the dynamical evolution is remarkable. The successful repro- 
duction of the analytical profile by the dynamically evolved 
particle distribution provides strong evidence for the capa- 
bility of the present model. There appear small deviations 
at the centre for the more concentrated hosts which may 
be caused by a lack of numerical resolution or by a slight 
inac curacy of the employed accretion formul a at early times 
(see lTasitsiomi et alj|200l ; IZhao et al.ll2003l ). However, this 
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Figure 3. Arbitrarily scaled number density profiles of the parti- 
cles for three host concentrations. The thin solid lines display the 
NFW-profiles for the corresponding concentrations. The group- 
centric distances are scaled by r s , the characteristic radius of the 
NFW-profile (r s = R vir /c). 



occurs only within a few per cent of the virial radius and 
will not affect the results presented here. 



3.3 Velocity dispersion profiles 

We now focus on the velocity profiles of the point mass 
sample. In particular we derive a general expression for the 
anisotropy parameter = 1 — (at,i/a r ) 2 where ov is the ra- 
dial velocity dispersion and <Jt,i is the uni-dimensional tan- 
gential velocity dispersion. If the tangential velocities are 
isotropic a t ,i = <r t /V2, where at is the two dimensional tan- 
gential velocity dispersion. By scaling the host-centric dis- 
tances by r s , the characteristic radius of the NFW-profiles, 
one achieves coincidence of the /3-profiles independent of the 
concentration. 

The upper panels in Fig. [4] display the radial mean ve- 
locities and the velocity dispersions of the particles for the 
three example hosts (c = 6, 10,20). Negative values in the 
fist panel indicate net infall. For the more concentrated hosts 
the radial velocities are close to zero, indicating fairly re- 
laxed systems. The particles within the cq host show some 
slight net infall at ~ 7r 3 which my be a result of the re- 
cent rapid growth phase of the host haloe. The velocity dis- 
persions in the secon d panel show a b ehaviour known from 
N-body simulations (|Hoeft et al.ll2004T ). they peak at about 
0.8r 3 . In the third panel the velocity dispersions are split up 
into radial and tangential components. To allow for a direct 
comparison with the radial component the uni-dimensional 
tangential velocity dispersion a t ,i is shown here. Finally, 
the anisotropy profiles /3(r/r s ) for the three host haloes are 
given in the lowest panel. Despite the different host con- 
centrations, we find that the anisotropy profiles are very 
similar if the distances scaled by the characteristic radius 
r s of the NFW-profile rather than the virial radius R v i r . 
The /3-profiles can be fit with an concentration-independent 
expression: 
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Figure 4. The upper three panels display mean radial velocity 
(v rad ), the velocity dispersion (a) and the orthogonal velocity 
dispersion components (radial oy and ID tangential cr t i) of the 
particles within host haloes of three concentrations. The group- 
centric distances are scaled by r s , the characteristic radius of the 
NFW-profiles (r a = R v i r /c). Negative values for the mean veloc- 
ities in the fist panel indicate net infall. The lowest panel shows 
the associated anisotropy profiles /3(r). The black solid line dis- 
plays the /3-fit given in Eq llll with an exponent of a = 0.45 and 
a scaling b = 0.55. Note, the anisotropy profiles of different con- 
centrations can be fitted equally well with the same parameters. 



This expression enables an analytical integration of the 
Jeans equations as discussed in the Appendix. 

In addition to the successful reproduction of the den- 
sity profiles the present inspection of the velocity profiles 
confirms that the dynamical properties of our model are 
also in good agreement with the dynamical characteristics 
of CDM haloes generated by N- body experim e nts. B ased 
on 10 relaxed dark matter halos IWoitak et al.l (|2005h find 
a similar shape and asymptotic values for the average /3- 
profile. Also the /3-profiles b ased on six highly resolved N - 
body systems presented by Dehncn & McLaughlin (2005) 



p(r) = b 



x + 1 



— r y r 



(11) 



(ba sed on simulations disc ussed in iDiemand et al.l (2004) 
and lDiemand et alj (l2004al )) show very similar features. Fi- 
nally, the density slope - velocity anisotropy relation dis- 
cussed in lHansen fc Moorel |2006) implies a self similarity 
of /3-profiles when scaled by r a as advocated here. 
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Figure 5. Mean velocity dispersions within the listed fractions 
of the virial radius versus the host concentration. The thick lines 
display the results from the dynamical models and the thin solid 
lines arise from the numerical integration of the Jeans equation 
using the /3-formula Eq. 1111 (see Appendix). 



3.4 Concentration-sigma relation 

We now investigate the relation between the mean particle 
velocity dispersion and the concentration of the host haloes. 
It is expected that the velocity dispersion increases with 
concentration, therefore the present results should be con- 
sidered as an further affirmation of the model and a point 
of reference for subsequent investigations of the subhaloes 
sample. 

Fig. [S]shows the dependence of the mean particle veloc- 
ity dispersions on the host concentrations. The mean values 
of the dispersions are computed within several fractions of 
the virial radius. The thin solid lines display the results from 
the analytical integration of the Jeans equation using Eg. 1111 
to describe the radial behaviour of the anisotropy parame- 
ter. Details of this integration can be found in the Appendix. 
The agreement between the model and the analytical com- 
putation is convincing. For very small concentrations < 5 
some deviations appear, but hosts haloes with these small 
concentrations are in the process of formation and the ap- 
plicability of the Jeans equation in this form may be ques- 
tioned. The ~ 5 per cent deviation for the mean dispersion 
within the full virial radius 1.00-R„i r probably results from 
the orderly accretion mode inherent in our model. Its traces 
are certainly not blurred out before the second turnaround 
passage (see also Fig. [1} . The main finding is that the veloc- 
ity dispersion for hosts with concentrations ~ 20 are a factor 
of 1.2—1.3 higher then the mean velocity dispersions of equal 
mass hosts with concentrations ~ 5. The mean central ve- 
locity dispersions increase somewhat stronger with the host 
concentration than those within larger radii. 

In summary, our comprehensive examinations of the 
point mass sample have shown satisfactory agreement with 
the results from N-body simulations, indicating that our 
simplified dynamical model can reproduce the main charac- 
teristics found by more elaborate N-body simulations. How- 
ever, we notice that the outskirts (> 0.6-R„i r ) may be slightly 
affected by the unphysical regularity of the accretion mode. 
Nevertheless, the model enables us to measure in detail the 



dependence of the subhaloe properties on the host concen- 
tration. 



4 RESULTS FOR THE SUBHALOE SAMPLES 

This section describes the dynamical evolution of extended 
subhaloes within the growing potentials of their hosts. As 
described in § 12.2.21 we attribute to every (sub)haloe a 
NFW-density profile and include the effects of dynamical 
friction and tidal truncation in the integration of the or- 
bits. For every host we distinguish two subhaloe samples 
at z — 0, a complete sample and a tidally reduced sample. 
The former comprises all 100000 subhaloes independently 
of their present truncation radius and the latter counts only 
those subhaloes whose tidal truncation radius is larger then 
3kpc. The tidally reduced sample is a sub-sample of the 
complete sample. For example, the reduced samples of the 
C6, Cio and C20 host haloes comprise a number of 74777, 
57218, 36258 surviving subhaloes, respectively, of the orig- 
inal 100000 within the complete sample. The idea behind 
this distinction is to imitate the effects of detection limits 
which are inherent to observation and simulations. If satel- 
lite galaxies fall below a certain luminosity, they won't ap- 
pear in the catalogue. And if subhaloes fall below a certain 
particle limit they won't be counted as substructure. In both 
cases the disappearance from the catalogue does not neces- 
sarily mean total physical disassembly. The accretion times 
and the assembly history of the subhaloes samples do not 
differ from those plotted in Figs. [1] and [2j for the point mass 
sample. Consequently, we start the investigation of the sub- 
haloe samples with the discussion of the evolution of the 
mass function. 

4.1 Mass functions 

The assignment of initial subhaloe masses is designed to 
reproduce a differential mass distribution dN/dM with a 
power law shape and a logarithmic slope of a — —1.86, the 
appropriate slope for field haloes measured in N-body simu- 
lations. The upper panel of Fig. [6] compares the differential 
mass functions of the arriving subhaloes (thin lines) with the 
mass function of the complete subhaloes at z = after they 
were dynamically processed within the three example hosts 
(c = 6, 10, 20). The mass functions are shown within a range 
of 5 x 10 s - 2 x 10 12 M . The initial mass functions for the 
three concentrations perfectly coincide, indicating that the 
temporarily slightly different upper mass limits due to the 
altering host growth rates (see Eq.[7]) do not affect the overall 
mass function. Interestingly the logarithmic slope does not 
change for the dynamically processed subhaloe samples at 
z = 0. Similar results have been found in N-body simulations 
(e.g. iReed et al.ll2005l ) . There appears a slight dependence of 
the amplitude of the dynamically processed mass functions 
on the host concentrations. More concentrated hosts have 
slightly reduced mass functions. The reason for the weak 
dependence on host concentration can be explained by the 
correlation between mean subhaloe accretion times and host 
concentrations as displayed in Fig. [7] On average subhaloes 
in more concentrated hosts are accreted earlier, thus they 
are exposed to the tidal field of the host for a longer time 
and they suffer more severe mass losses. The thick lines in 
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Figure 6. The upper panel compares the differential mass func- 
tion of the subhaloes in the complete sample before entering the 
virial radius of the host (thin lines) with the mass function of the 
subhaloes at z = (thick lines). The straight line displays the 
slope used in Eq. [7] There is a week dependence of the amplitude 
of the final mass function on the concentration. The lower panel 
shows the ratios of the initial to the present mass functions. Thick 
lines are used for the complete subhaloe sample (these values can 
be directly read off the upper panel). The thin lines in the lower 
panel display the ratios between the mass function of the initial 
subhaloe mass function to the tidally reduced sample. 



the lower panel of Fig.[6]may clarify this issue. They display 
the ratios of initial to the final subhaloe mass function. The 
C6 host exhibits a difference in amplitude of a factor of ~ 2, 
whereas in the C20 host the final mass function is reduced by 
a factor of ~ 3. The thin lines in the lower panel display the 
ratios between the initial and the final mass distributions 
in the tidally reduced subhaloe sample. At the low mass 
end the difference between initial and final mass function is 
enhanced, however the high mass end shows barely any dif- 
ference. Therefore, at the low mass end the power law slope 
is conserved even in the tidally reduced sample. The solid 
line in Fig. [7]displays the mean accretion times as a function 
of host concentration for the tidally reduced sample, fn par- 
ticular for concentrations ~ 20 the difference of the mean 
accretion time is about 3.5 Gyr. 



4.2 Density profiles 

Fig. H] shows the number density profiles for the complete 
and the tidally reduced subhaloe samples within the three 
host haloes with present concentrations of 6, 10 and 20. In 
contrast to the density profiles in Fig. [3] the distances from 
the centre of the host are scaled by the virial radius i?„ir- 
For comparison the associated NFW-profiles are displayed 
with thin solid lines. Similar to findings for the point mass 
sample in § 13.21 the agreement between the NFW-profile and 
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Figure 7. Variation of mean accretion redshift of subhaloes with 
host concentration at z = 0. The dotted line shows the results for 
the complete subhaloe sample. The solid line displays the relation 
for the tidally reduced subhaloe sample. 
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Figure 8. Arbitrarily scaled number density profiles for the com- 
plete and the tidally reduced subhaloe sample within host haloes 
of three concentrations (c=6,10,20). The thin solid lines give 
the associated NFW-profiles. The complete sample shows per- 
fect agreement with the NFW-profiles. The profiles of the tidally 
reduced samples flatten towards the centre. Subhaloe populations 
in more concentrated hosts are more efficiently reduced. 



the complete sample is remarkable, indicating that the com- 
bined action of dynamical friction and tidal truncation does 
not severely alter the number density profile. Subhaloes lose 
a substantial fraction of their mass during their first peri- 
centre passage, after that the dynamical friction is strongly 
reduced since dynamical friction drag is proportional to the 
subhaloe mass (Eq. 1 10 p . The number densities of all tidally 
reduced samples show a centrally flattened behaviour, in- 
deed at the very centre they display positive gradients. We 
do not discuss to what extent these positive gradients are 
physical since the presence of baryons in subhaloes certainly 
alters their tidal resistance. However, it is noteworthy that 
most of tidal reduction within the ce host is confined to the 
volume within half the virial radius whereas the tidal reduc- 
tion is visible over the whole haloe volume in the case of the 
C20 host. 
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Figure 9. The upper three panels display mean radial velocity 
( v rad)} the velocity dispersion (a) and the orthogonal velocity 
dispersion components (radial o> and ID tangential 0^,1) of the 
complete subhaloe population within host haloes of the indicated 
concentrations. As distinct from Fig. [4] here the group-centric 
distances are scaled by virial radius R v i r . 
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Figure 10. The upper three panels display mean radial veloc- 
ity (v r ad)i * ne velocity dispersion (<r) and the orthogonal velocity 
dispersion components (radial o> and ID tangential cr ti i) of the 
tidally reduced subhaloe population within host haloes of the in- 
dicated concentrations. As distinct from Fig. [4] here the group- 
centric distances are scaled by virial radius R v i r . 



4.3 Velocity dispersion profiles 

Some very interesting insights into the dynamics of sub- 
haloes can be gained by comparisons of the velocity dis- 
tributions for the complete and the tidally reduced samples 
displayed in Figs. l9l and 1101 respectively. It is important to 
keep in mind that the tidally reduced sample is a subsample 
of the complete sample. Therefore, higher mean velocities 
in the reduced sample can only be achieved if preferentially 
slow moving haloes are removed. In other words there is no 
source of high velocity haloes in the reduced sample. The 
uppermost panel in Fig. [9] displays the mean radial velocity 
of the complete sample. With exception of the outskirts they 
are close to zero, indicating a relaxed state. On the other side 
the mean velocities of the tidally reduced samples (Fig. I10[) 
show a clear trend for inward motion, negative v ra d- This 
is most naturally explained by tidal shrinkage as subhaloes 
first approach the host centre. Consequently, most probably 
the truncation radii falls below the rejection radius of 3kpc 
close to the peri-centre passage of the subhaloe, resulting in a 
lack of outward streaming subhaloes in the reduced sample. 
For large radii the velocity dispersions of the two samples 



are similar (second panels in Figs. [9] and 1 10 p . However, the 
central velocity dispersions of the reduced sample are ~ 1.3 
times larger than in the complete sample. Again this increase 
of the velocity dispersion can only be explained by the lack 
of preferentially slow moving subhaloes. The comparison of 
the orthogonal components of the velocity dispersions (third 
panel in Figs. [9] and I10p indicates that the increase of the 
central velocity dispersion in the reduced sample is based on 
the strong enhancement of the tangential dispersions, which 
again must arise due to a lack of slow tangential motions. 
The radial components give the reverse picture. The central 
radial velocity dispersion of the reduced sample is actually 
smaller, indicating a deficit of fast radially moving subhaloes 
compared to the complete sample. A hint f or this mechanism 
was reported bv lFaltenbacher et al.l (|2005h investigating the 
velocity distribution of galaxies in hydro-dynamical simula- 
tions of clusters of galaxies. The resulting anisotropy profiles 
are quite different (note the different scaling in the lowest 
panels of Figs.l^land llOp . If one assumes similar processes op- 
erating for satellite galaxy population s in simulated groups , 
the negative values for j3 displayed in lBenatov et al] (|2006f ) 
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Figure 11. Variation of mean subhaloe velocity dispersions with 
concentration of the host for various inclusion radii as listed. The 
upper line of each pair of lines with the same line style always 
displays the dispersion of the tidally reduced sample. The lower 
lines show the values belonging to the complete sample. 

(thei r Fig. 5) are in agree ment with the findings presented 
here. iMaccio et al.l (|2006h discussed the enhanced survival 
rate of subhaloes with condensed baryonic cores. In theory 
the behaviour of the velocity distribution of satellite galax- 
ies in groups can be used to determine their tidal disruption 
rate. 

4.4 Concentration-sigma relation 

In analogy to FigO we display in Fig llll the mean disper- 
sions within 0.2, 0.4, 0.6 and 1 times the virial radius R v i r 
for the complete and the tidally reduced subhaloe samples. 
The larger dispersions always correspond to the tidally re- 
duced sample whereas the lower line of each pair refers to the 
complete sample. The dispersions of the point mass sample 
in Fig[S]are quite similar to the values of the complete sam- 
ple, indicating only a weak impact of dynamical friction on 
the global properties of the subhaloe population. However, 
the impact of tidal reduction is enormous for the central 
dispersions, e.g. for host concentrations of ~ 20 the veloc- 
ity dispersions within 0.2_R„i r are enhanced by a factor of 
1.3 compared to the complete sample. The dependence on 
concentration within the tidally reduced sample is also more 
prominent, the central dispersion in hosts with concentra- 
tion ~ 5 to hosts with concentration ~ 20 increases by a 
factor 1.4. The mean dispersion within the virial radius is 
only weakly affected by the tidal reduction. 

Since the complete sample resembles the behaviour of 
the point mass sample, the Jeans equation holds also in this 
case. Due to the negative mean radial velocities in the re- 
duced sample (Fig. Hip one of the preconditions for the de- 
scription by the Jeans equation is no longer valid (see Ap- 
pendix) . We do not aim to present a modified solution which 
can cope with radial motions, because the extent of inflow 
depends on the choice of the rejection radius in our model. 

The important result here is that an increase in cen- 
tral velocity dispersion of satellite galaxies in groups caused 
by higher concentration may be amplified if some fraction 
of the satellite population has been tidally disrupted. Both 



trends are in the same direction. Low concentration hosts 
have relatively low central velocity dispersions and their re- 
cent formation times make tidal disruption as a driver for 
an enhanced velocity dispersions unlikely. Strongly concen- 
trated hosts have higher velocity dispersions and accreted 
most of their satellites early on. This old population may be 
more severely reduced by tidal forces, causing an additional 
boost to the central velocity dispersion. 



5 SUMMARY 

Using a dynamical model for the evolution of subhaloes 
within a growing host potential in galaxy groups, we have in- 
vestigated the relation between host concentration and prop- 
erties of subhaloe population at * = 0. The concentration is 
correlated with the group formation redshift. Groups with 
equal virial masses at z — but different formation times 
show significant differences in their subhaloe properties. As 
a base point we have computed the dynamical evolution of 
point mass samples and their dependence on the host con- 
centration. The model then has been adapted to cope with 
the spatial extent of subhaloes and the impact of dynami- 
cal friction and tidal disruption. We only focus on the dy- 
namical evolution of subhaloes that are massive enough to 
possibly host galaxies. Every substructure analysis is sub- 
ject to detection limits whether it is the sensitivity limit in 
connection with observations or the resolution limit in nu- 
merical simulations. To study the effects of detection limits, 
we create a reduced subhaloe sample comprising only those 
subhaloes with current radii larger than a rejection radius 
of r re j = 3kpc. The reduced sample is a sub-sample of the 
complete sample, comprising all subhaloes independent of 
the actual truncation radius. The rejection radius can be 
considered as a parameter of the model and has been cali- 
brated by observations as well as numerical simulations (see 
iFaltenbacher fc Mathewsll2005t ). Originally the model is set 
up to trace the evolution of the dark matter subhaloes in 
a group environment, however if the subhaloes host lumi- 
nous galaxies, a tidal truncation below 3 kpc most probably 
would affect the stellar component of the galaxy as well. 

Our main findings are: (1) The velocity dispersion of 
subhaloes or satellite galaxies in groups depends on the 
concentration of the underlying dark matter distribution of 
the host haloe. Equal mass hosts with higher concentra- 
tions have enhanced velocity dispersions of the subhaloes, 
in particular at the centre. Since higher host concentrations 
are caused by earlier formation times, this relation implies 
that subhaloe populations with early formed hosts exhibit 
higher velocity dispersions than those residing in more re- 
cently formed hosts of equal mass. (2) We propose a fitting 
formula for the radial dependence of the anisotropy param- 
eter j3, which is independent of concentration if the group- 
centric distances when scaled by the characteristic radius of 
the NFW-profile r s . (3) The increase of the central veloc- 
ity dispersion with concentration is amplified in the reduced 
sample. This is caused by an enhancement of the tangential 
velocity dispersion and a lack of slow radial motions com- 
pared to the complete sample. (4) The dynamical evolution 
of subhaloes that experience tidal stripping and dynami- 
cal friction does not alter the slope of the mass function, 
in agreement with the results from N-body simulations (see 
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e.g. lReed et~alll2005h . but dynamical evolution reduces the 
amplitude by a factor of ~ 2. In the complete sample the 
amplitude of the mass function is only weakly dependent 
upon the host concentration. However, the reduced sample 
shows a stronger dependence on host concentration with am- 
plitudes decreasing with increasing host concentrations. For 
a host concentration of 20 the amplitude of the initial mass 
function decreases by a factor of ~ 5, however this value 
depends strongly on the choice of the rejection radius. 

The present investigation reveals the difficulties asso- 
ciated with mass estimates derived from the velocity dis- 
persion of satellites in groups. The concentration of the 
host haloe and the effects of tidal reduction of the satel- 
lite galaxies in groups may alter the central velocity disper- 
sions. The appearance of intra-group stars is strong evidence 
for the reduction of the sa tellite luminosities. According to 
lOsmond fc Ponmanl l|2004h the velocity dispersion provides 
a very unreliable measure of system mass. The concentra- 
tion and its influence on tidal reduction add more concerns 
on the accuracy of such mass estimates. 

Combined observations of X-ray temperatures and ve- 
locity dispersions of satellite galaxies have recently become 
available for groups. As a application of the present analysis 
these observations my be used to infer the impact of tidal 
forces onto the central satellite population. If tidal forces 
in fossil groups efficiently reduce the number of detectable 
satellites, we expect the central velocity dispersion to be 
high compared to the X-ray te mperature. Some support f or 
this scenario has been found bv lKhosroshahi et all [2006 V 

A comparison of the X-ray temperatures of the hot 
intra-group gas with the velocity dispersions found for the 
satellite galaxies can be used to infer the degree of tidal dis- 
ruption. A flattening of the central number density profiles 
along with an increasing velocity dispersion are the signa- 
tures of tidal reduction of the central subhaloe or satellite 
pop ulation. 

ISommer-Larsenl (|2006l ) considers radial velocity disper- 
sions in fossil systems at a fixed radius r > 30 kpc. We have 
presented a uniform fitting formula for the anisotropy pa- 
rameter provided the group-centric distances are scaled by 
r s , the characteristic radius of the NFW-profile of the host. 
Fos sil groups are assume d to be early formed systems (see 
e.g. iD'Onghia et"ai1l2005h with high concentrations. Our fit- 
ting formula implies that at a fixed radius more concentrated 
hosts generate more radially anisotropic velocity dispersions. 
This result has to be taken into account if the velocity dis- 
persions of different groups are compared at physically sim- 
ilar radii. 
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For completeness we repeat the fitting formula for the 
anisotropy parameter f3 as given in Eq. [11] 



I3{r) = b 



x + 1 



x = r/r s , (15) 



where r a is the scale factor of the NFW-profile. With this 
expression for /3(r) the integrating factor (Eq. I13[) can be 
written as 



u(x) = u(xr s ) — exp {26 2-F1 [a, a, 1 + a, —x]} 



(16) 



where 2F1 is the hyper-geometric function and x = r/r s 
as introduced in the foregoing equation. Finally, the mass 
weighted mean velocity dispersion within the radius r is 
given as 



a{r) = 



£ p{r')al{r')^dr' 
fa p(r')r ,2 dr' . 



(17) 



The mean velocity dispersions represented by the thin lines 
in Fig. [5] are computed on the assumption that the matter 
distribution follows a NFW-density profile with a virial mass 
of 3.9 x 10 13 M . As seen in Fig. □ up to 10 per cent of all 
the particles are located outside the virial radius at z = 0. 
This particle deficit at large radii may explain the ~ 5 per 
cent deviations between the velocity dispersions obtained 
from the model and the Jeans equation (r < 1.00R v i r line 
in Fig. [5j) . 



APPENDIX: INTEGRATION OF THE JEANS 
EQUATION 

Fig. [S] compares the velocity dispersions from our dynam- 
ical model with the predictions of the Jeans equation. For 
a spherical symmetric and static (oy = at = 0, no mean 
radial or tangential motions) system the Jeans equation can 
be written as 



~-rG°°v) + 2/3 — = — j- 
p dr r dr 



(12) 



where p is the density, oy is the radial velocity dispersion, $ 
is the gravita tional potential and P is th e anisotropy param- 
eter (see e.g. iBinnev fc Tremaine! [l987). A solution of this 
first order differential equation can be be obtained by means 
of multiplying both sides of the Eq. [12] by the integrating 
factor 



u(r) 



exp 



r /3(r')dr' 



(13) 



and performing an integration from a given radius r to infin- 
ity. The application of the boundary condition lim, — >00 ay = 
which any ordinary bound system satisfies results in 



p(r)<j r {r) exp 



r /3(r')dr' 



f 



— (r )p{r ) exp 



P{r")dr" 



dr' . (14) 



